Controlling surface morphologies by time-delayed feedback 
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We propose a new method to control the roughness of a growing surface, via a time-delayed 
feedback scheme. As an illustration, we apply this method to the Kardar-Parisi-Zhang equation in 
1+1 dimensions and show that the effective growth exponent of the surface width can be stabilized 
at any desired value in the interval [0.25,0.33], for a significant length of time. The method is 
quite general and can be applied to a wide range of growth phenomena. A possible experimental 
realization is suggested. 
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Introduction. The control of unstable states in chaotic 
or pattern-forming nonlinear dynamic systems has at- 
tracted much interest recently [H, 0] . Time-delayed feed- 
back control Q has been especially successful in stabiliz- 
ing a variety of dynamic and spatial structures, includ- 
ing noise-induced oscillations andpatterns found, e.g., 
in semiconductor nanostructures [J H, H, 0] • The fabri- 
cation of such nanostructures typically involves the de- 
position of a material onto a substrate. One of the pri- 
mary experimental goals is to achieve nanoscale control 
of layer thickness and surface (or interface) morphology. 
On the theoretical side, considerable effort has focused 
on developing suitable evolution equations for the grow- 
ing layer and its surface [g| . While many different ver- 
sions % M, HO, E 03 M, EE El of these equations 
exist, depending on the details of deposition processes 
and molecular interactions and kinetics, all of them share 
certain fundamental characteristics: they are noisy, non- 
linear partial differential equations in space and time, 
and describe an important class of generic nonequilib- 
rium phenomena. 

It is natural to ask whether the control techniques of 
nonlinear dynamics can be successfully applied to sur- 
face growth problems. The goal is, of course, to sta- 
bilize desired surface characteristics, such as its spatio- 
temporal height-height correlations or its roughness dur- 
ing the growth process. Even if such control can only 
be achieved in a finite window of time, its experimen- 
tal potential is undiminished since the deposition process 
can simply be terminated at the desired time, thanks to 
today's precise in situ characterization capabilities. In 
this letter, we provide a first set of answers to this ques- 
tion. We choose the most promising type of control, time- 
delayed feedback, and study its effects on a paradigmatic 
growth model, the Kardar-Parisi-Zhang (KPZ) equation 
[17( . Specifically, we attempt to control the effective dy- 
namic growth exponent (3 associated with the roughness 



of the growing surface. Implementing two realizations of 
the control scheme, we will see below that we can indeed 
stabilize (3 in a range of values between the two universal 
limits, 1/4 and 1/3, over at least one to two decades in 
time. 

This letter is organized as follows. We first introduce the 
KPZ equation and our numerical solution scheme, ac- 
companied by a representative data set for the growing 
surface roughness, and specifically its growth exponent, 
in the absence of control. Next, we implement two types 
of time-delayed feedback control and demonstrate how 
reasonably accurate values of the (effective) growth ex- 
ponent can be achieved. We conclude with a summary 
and a discussion of open questions. 
The KPZ equation. If overhangs and bulk fluctuations 
can be neglected, growth phenomena are often modelled 
in terms of nonlinear stochastic partial differential equa- 
tions. A single- valued variable, h(x, t), denotes the height 
of the surface above a reference plane and fluctuates as 
a function of time t and position x (measured in this d- 
dimensional plane). The simplest such equation is the 
KPZ equation [ItJ which describes the growth of a sur- 
face in the absence of any conservation laws: 

d t h{x,t) = vV 2 h{x,t) + ^{Vhf + V (x 1 t) (1) 

Here, v > denotes an interface smoothing term, asso- 
ciated with a surface tension; the nonlinear coupling A 
reflects the strength of lateral growth, and rj(x, t) mod- 
els the height fluctuations due to random deposition of 
material. An overall, constant growth velocity has al- 
ready been eliminated, by transforming into a suitable co- 
moving frame. Focusing on large-scale, long-time prop- 
erties of the surface, it is sufficient to consider Gaussian 
white noise, i.e., 



< r](x,t)ri(x',t r ) > 





2DS d (x - x')S(t - t') (2) 
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The KPZ equation has been discussed in man y d iffer- 

ent contexts, including thin film growth 

fluctuating hydrodynamics [22j, driven diffusive systems 
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|23l.l24j]. tumor growth in biophysics 1251, 12611 . propagating 
fire fronts [2?], and econophvsics [29j| In the follow- 
ing, we will use the language of surface growth, but our 
findings are easily translated into these other contexts 
and just as relevant there. 

For simplicity, we restrict ourselves to one spatial dimen- 
sion. We monitor the time-dependence of the (root mean 
square) surface roughness w, defined by 



(3) 



Here, L denotes the system size, and h(t) = 
L^ 1 h(x, t) is the mean surface height at time t. Con- 
figurational averages are denoted by (...). The sum over 
x anticipates the space discretization associated with the 
numerical integration scheme. It is well known that w 
obeys scaling in the form w(L,t) = L a f(t/L z ) where / 
is a scaling function and a and z denote the roughness 
and dynamic exponents, respectively [30j] . In the satura- 
tion regime t/L z 3> 1, / approaches a constant so that 
w ~ L a becomes independent of time. In contrast, in 
the growth regime t/L z <§; 1, the width grows as a power 
of time, w ~ t@ , with a growth exponent (3. Consistency 
with the general scaling form imposes an exponent iden- 
tity, (3 — a/z. The values for the exponents a, (3, and 
z are universal. Two universality classes can be distin- 
guished. If the nonlinear term of the KPZ equation van- 
ishes (A = 0), the equation reduces to the exactly soluble 
Edwards- Wilkinson (EW) equation [3l[, with a = 1/2 
and z = 2 whence (3 = 1/4. In contrast, any nonzero 
value of A belongs to the KPZ universality class with 
a = 1/2 and z = 3/2 whence j3 = 1/3 [13]. These val- 
ues give us some benchmarks against which we can check 
our numerical scheme. We use a forward-backward Eu- 
ler method [HJ to solve the stochastic differential equa- 
tions numerically. The associated discretization captures 
certain prefactors much more accurately than the stan- 
dard (naive) scheme 33] . Our integration parameters are 
v = 0.1 and D = 0.5; the space discretization is set to 
Ax = 1, with L = 1024 and 4096, and the time increment 
is set at At — 10~ 3 . Finally, A varies between 0.00 and 
0.25. The upper cutoff is chosen so as to avoid numerical 
instabilities. 

Figure [1] shows a scaling plot for w(L,t) before any con- 
trol schemes are implemented. Data for three differ- 
ent values of A are shown. The roughness exponent a 
is consistent with 1/2, independent of A, as expected. 
For A = 0, we see excellent data collapse with the EW 
scaling exponents, and the KPZ exponents are confirmed 
for the largest A. The latter should be universal, for all 
A ^ 0; however, for < A < 0.25, strong crossover effects 
between EW and KPZ behavior are observed. Remark- 
ably, this crossover manifests itself as a surprisingly clean 
power law, with a A-dependent effective growth exponent 
(3 below 1/3. Eventually, the asymptotic value (1/3) is 
reached, but only after an L- and A-dependent crossover 
time. 
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Figure 1: (color online). Scaling plots of the rough- 
ness w(L,t) and height-height correlation function C(l) = 
(h(l,i)h(0,i)} ~ l a (insets) for system sizes L = 1024 (red) 
and L = 4096 (black) and three choices of A: (i) A = (shifted 
by a factor 100), (ii) A = 0.1 (shifted by a factor 10), (iii) 
A = 0.25. v = 0.1 for all data sets. The green broken lines 
provide guides to the eye. 



Systems with control. We now turn to possible con- 
trol mechanisms. For chemical vapor deposition of silica 
films, there is some experimental evidence [l8[ that the 
lateral growth velocity is related to the temperature, via 
a temperature-depending sticking probability. In other 
words, A - and hence effective growth exponents - can 
be controlled via the temperature. For our differential 
equation, we tune A directly, in order to stabilize a de- 
sired effective growth exponent, (3q. In detail, the scheme 
is as follows. First, we choose the desired value of the 
growth exponent, /3o, and select an appropriate time de- 
lay r. Generating sufficiently many samples of h(x, t), we 
record w(t — r) and w(t) (the argument L will be omit- 
ted from now on) . The local exponent (3i oca i at time t is 
defined as 



(3local{t) 



logw(i) — \ogw{t — t) 
\ogt - \og{t - t) 



(4) 



Depending on the sign and value of (3i oca i{t) — (3q, we 
adjust the nonlinear coupling, A, of the KPZ equation, 
as follows. First, we introduce a control function F(t). 
For digital control, we define 



F(t) 



if Plocal < (3q 
if Plocal > (3o 



(5) 



where the parameter a defines the control "bit", i.e., the 
amount by which A changes at each control step. Al- 
ternatively, we also investigate a differential method for 



which 

F(i) EE K(fio - Plocal) (6) 

and K sets the amplitude of the control strength. Given 
one of the two choices of F(t), the control scheme sets in 
at time to- From then on, the nonlinearity A is updated 
at times t n = t + nr, n = 1,2,..., starting from an 
initial value Ao , according to 

{A , ift<t 
X(t-T)+F(t), ift = t n (7) 
A (i n )) if i„ < i < t n+ i 

Our scheme is successful if Piocai(t) approaches (3q and 
then settles at the desired value within a reasonable time 
frame after the control has been activated. 
Some comments are in order. Starting from a random 
initial condition, we first choose a starting value, Ao, for 
the nonlinearity and integrate the KPZ equation without 
control up to time to, in order to eliminate transients. 
A reasonably stable growth regime is achieved around 
to ~ 10, independent of L (provided L is not too small, 
i.e., L > 64). Then, we turn on the control, following 
either the digital or the differential scenario. Regarding 
the choice of the time delay r, it must be large enough 
compared to the time increment At so as not to interfere 
with the integration procedure, but small enough to pro- 
vide responsive control. We find that we get good results 
for a time delay 0.1 < t < 1.0. Similarly, we choose the 
control amplitudes a and K such that the increments in A 
are small compared to Ao but large enough to generate a 
noticeable response. For example for r = 1.0, choosing a 
in the range [0.002,0.02] and K in the range [0.005,0.05] 
provides the best results. 

Results. Figures [2] and [3] show our results. Starting 
from three initial values of Ao, namely, 0, 0.1, and 0.25, 
we attempt to stabilize the effective growth exponent at 
00 = 0.29, mid-way between the KPZ and EW values. Ir- 
respective of Ao, we find that both digital and differential 
control result in an effective growth exponent very close 
to the desired value, over at least a decade of integration 
time (100 < t < 1000). 

For sufficiently large time, the control function A(t) ap- 
pears to approach a constant value, close to 0.15. How- 
ever, the details of the approach depend on the initial Ao- 
For strong initial nonlinearity Ao = 0.25, X(t) is approx- 
imately monotonically decreasing, apart from significant 
fluctuations. For both weak and vanishing initial cou- 
pling, X(t) approaches its "limit" from below. This be- 
havior is observed for both, digital and differential, con- 
trol. We tested several other choices of 0.25 < flo < 0.33 
and Ao, and found similar behavior. 
Experimentally, it is usually desirable to achieve small 
roughness. To push our control schemes to the limit, we 
test the most extreme case, namely, Po = 0.25 with large 
initial nonlinearity Ao = 0.25. Before the control sets in, 
the roughness grows considerably faster than t ' 25 . As 
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Figure 2: (color online). Roughness and control function (in- 
set) evolution for digital time delayed feedback control for 
strong (blue), weak (red), and zero (black) initial nonlinearity 
Ao. The desired effective growth exponent is set at /3o = 0.29. 
To provide a comparison, the straight (green) lines have slopes 
0.33 (dotted), 0.29 (dashed), and 0.25 (dash-dotted). All data 
sets are obtained with v — 0.1 and a — 0.01. 
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Figure 3: (color online). Roughness and control function (in- 
set) evolution for differential time delayed feedback control for 
strong (blue), weak (red), and zero (black) initial nonlinearity 
Ao. The desired effective growth exponent is set at /3o = 0.29. 
To provide a comparison, the straight (green) lines have slopes 
0.33 (dotted), 0.29 (dashed), and 0.25 (dash-dotted). All data 
sets are obtained with v — 0.1 and K — 0.02. 



soon as the control sets in, X(t) decreases quite dramati- 
cally, leading to a reduction of the effective j3. However, 
over the time period considered (i < 1000), it never de- 
creases far enough to reach the desired 0.25. 
Finally, we note that it is not possible to achieve ex- 
ponent values outside the interval [0.25,0.33]. Choosing 
/3 > 0.33 generates unbounded growth of the control 
function A(i), accompanied by instabilities in the integra- 
tion routine. Similarly, (3o < 0.25 quickly leads to large 
negative values of X(t) which tend to favor KPZ expo- 
nents (since the sign of A plays no role). As a result, X(t) 
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becomes even more negative until a numerical instability 
occurs. To avoid this instability, we also implemented 
a symmetrized version of control (with a —> —a when 
A(i) < 0). In this case, X(t) approaches zero and fluc- 
tuates about it, so that the effective exponent settles at 
0.25. 

Conclusions. To summarize, both digital and differen- 
tial control are rather successful at stabilizing effective 
growth exponents in the KPZ equation. For the rela- 
tively small system sizes used here, these exponents can 
be tuned in the range [0.25,0.33], i.e., within the limits 
set by the EW and the KPZ equation, respectively. Let 
us emphasize again that only the values 1/4 (for A = 0) 
and 1/3 (for any A / 0) correspond to true asymptotic 
exponents; for larger system sizes and longer integration 
times, these emerge clearly. However, for small systems, 
we observe surprisingly clean effective growth exponents 
which appear to depend monotonically on the magni- 
tude of the nonlinearity. Hence, it is possible to choose 
a desired reference exponent [3q and implement a time- 
delayed control of the nonlinearity in such a way that the 
effective exponent first approaches (3q and then stabilizes 
at that value for a significant length of time (roughly, 
10 2 ^ t < 10 3 , in our units). In all simulations, the 



(stationary) roughness remains constant at a = 0.5, re- 
flecting the value 1/2 which is common to both the EW 
and the KPZ equation. The control protocol itself is in- 
dependent from the dimension of the surface. Work is in 
progress to test other growth equations and to extend the 
KPZ study to 2+1 dimensions. We thus hope to develop 
a control tool which might also be useful in experimen- 
tal setups. The work of Ojeda et al [H, [l9| gives some 
indications that the KPZ nonlinearity is tunable via the 
temperature. Hence, it seems feasable to implement a 
time-delayed feedback loop and stabilize desired growth 
exponents by suitable adjustments of the temperature. 
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